Turbulent Vortex Flow Responses at the AB Interface in Rotating Superfluid 3 He-B 
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In a rotating two-phase sample of 3 He-B and magnetic-field stabilized 3 He-A the large difference in 
mutual friction dissipation at 0.20 T c gives rise to unusual vortex flow responses. We use noninvasive 
NMR techniques to monitor spin down and spin up of the B-phase superfluid component to a sudden 
change in the rotation velocity. Compared to measurements at low field with no A-phase, where 
these responses are laminar in cylindrically symmetric flow, spin down with vortices extending across 
the AB interface is found to be faster, indicating enhanced dissipation from turbulence. Spin up 
in turn is slower, owing to rapid annihilation of remanent vortices before the rotation increase. As 
confirmed by both our NMR signal analysis and vortex filament calculations, these observations are 
explained by the additional force acting on the B-phase vortex ends at the AB interface. 

PACS numbers: 67.30. he, 67.25.dk, 47.37.+q, 03.75.Kk 



Introduction:-A two-phase sample of superfluid 3 Hc 
is a unique case of coherent quantum matter, consisting 
of coexisting regions of magnetic-field stabilized 3 Hc-A 
and of 3 He-B at zero or low field. Here two coherent 
states of the same orbital (L = 1) and spin (5 = 1) 
triplet order parameter manifold can be investigated in 
phase equilibrium. The structure and dynamics of vor- 
tices at the AB interface have been studied earlier in a 
long rotating cylinder We use this setup to examine 
the responses of B-phase vortices to a step change in the 
rotation velocity Q. At low temperatures existing vor- 
tices become easily unstable in rapid changes of rotation 
0. As a result, both the formation of new vortices and 
the annihilation of existing vortices are associated with 
low energy barriers. A sudden increase (reduction) in 
leads to a continuous, but slow spin up (spin down) in the 
rotation of the superfluid component. In the two-phase 
sample at 0.20 T c mutual friction dissipation is estimated 
to be two orders of magnitude larger in the A phase [1] . 
As a result, A-phase vortex flow is almost instantaneous 
while the B phase responds only slowly. This leads to un- 
usual dynamics which differs remarkably from the lami- 
nar vortex flow observed in the absence of the A-phase 
stabilization field. 

In a cylinder filled with superfluid 4 He, spin up and 
spin down are assumed to be turbulent at essentially all 
temperatures below T\. In 3 He-B the responses have 
been found to be laminar for axially homogeneous spin 
down or spin up at least to below 0.20 T c However, 
in the two-phase sample B-phase spin down is faster and 
axially inhomogencous. NMR line shapes, which deviate 
from those measured for laminar vortex flow, can here 
be studied in controlled conditions in the cylindrically 
symmetric "flare-out" order parameter texture; in other 
words it is not broken cylindrical rotation symmetry, but 
the changed boundary conditions that cause the faster 
response and the increased dissipation. Thus the possi- 
bility to stabilize the A-phase layer is similar to changing 



in situ the boundary conditions at one of the end plates 
of the rotating cylinder. The resulting changes in the dy- 
namics indicate that boundary conditions arc important 
and provide further evidence for the fact that pinning at 
the end plates is responsible for the turbulent responses 
of superfluid 4 He. The 4 He vortex core diameter is of 
atomic size and is apparently strongly pinned at most 
surfaces, while in 3 He superfluids the core diameter is at 
least two orders of magnitude larger. 

Experimental techniques: -The two-phase liquid 
3 He sample is contained in a smooth-walled quartz cylin- 
der which is mounted on a nuclear demagnetization cool- 
ing stage and can be rotated at an angular velocity fi, by 
rotating the cryostat (Fig. []}. The time evolution and 
the spatial distribution of vortices is surveyed with two 
NMR detector coils. Two quartz tuning fork oscillators 
are included for temperature measurement. A small su- 
perconducting solenoid around the cylinder provides the 
axially oriented magnetic field for stabilizing an A-phase 
layer which divides the NMR sample in two identical B- 
phasc sections. 

Vortex structures in the A and B phases are incom- 
mensurable, but they do interconnect across the AB in- 
terface after the doubly-quantized A-phase vortices disso- 
ciate and terminate in point defects on the phase bound- 
ary The critical velocity of vortex formation is low 
in the A-phase section Q and in rotation it is approxi- 
mately in the equilibrium vortex state, where both the 
normal and superfluid components are in solid-body ro- 
tation with the container, v n = fi x r ~ v s . In the 
B-phase sections the critical velocity is an order of mag- 
nitude higher Q . Established procedures exist for main- 
taining in the B-phase sections either vortex-free counter- 
flow, where the superfluid fraction is at rest in the labora- 
tory frame (v s = 0), or a flow state with a central vortex 
cluster with a known number of vortex lines N. The clus- 
ter is formed from rectilinear vortices at the solid-body- 
rotation density 2S1/k, so that within the cluster v s m v n 
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FIG. 1: Measuring setup. Two independent cw NMR spec- 
trometers monitor vortex motion in the 110 mm long section 
of the sample cylinder above the 0.75 mm orifice. The mid- 
dle section with the quartz tuning fork oscillators is employed 
for thermometry. The bottom section below the 0.3 mm 
orifice provides thermal contact to the refrigerator. The su- 
perconducting solenoid in the center of the NMR sample pro- 
duces the axially oriented barrier field Hb for stabilizing 3 He- 
A. With a current of 6.2 A in the barrier magnet the width of 
the A-phase layer is 10 mm. 



while outside v s = k N / (2irr) (k = H/ (2777,3) is the 3 He 
circulation quantum). Thus there exist different configu- 
rations from where spin-down/spin- up measurements can 
be started. These can be classified according to how the 
51 range of the measurement interacts with the AB in- 
terface instability, which usually is characterized by its 
critical angular velocity 51ab (T, P) measured in a situa- 
tion when the A phase is in the equilibrium vortex state 
and the B phase is vortex- free [§] . 

Our spin-down measurement starts from the equilib- 
rium vortex state (at f2;) which extends through all three 
sample sections, with vortices interconnected across the 
two AB interfaces. Choosing Oi < Hab(T, P), the AB 
interface instability plays no role in these results. With 
decreasing temperature the 3 He-B mutual friction dissi- 
pation a(T) drops increasingly below that of 3 He-A and 
its characteristic dynamic response time, t(T) cx l/a(T), 
is slowed down. The faster A-phase dynamics provides 
then additional pull on the B-phase vortex ends on the 
phase boundary. This leads to faster and axially inhomo- 
geneous B-phase spin down (Fig. [2]). The ensuing events 
are the following: (i) The A-phase vortices spiral rapidly 
in laminar motion to the cylindrical wall. This exerts a 
force on the B-phase vortices such that the vortex den- 
sity in the center of the cylinder is depleted faster than at 
larger radii r and at larger distances Az from the AB in- 
terface. The central depletion in vortex density proceeds 
to such extent that it can be characterized as a vortex- 
free dome on the AB interface. Thus the vortex density 



FIG. 2: Calculated vortex configuration during spin down of 
the superfluid component, 30 s after a step-like reduction of 
Q, from 0.25rad/s to zero. (Left) View from the top into the 
cylinder showing the empty hollow in the center and twisted 
vortices around it. (Right) Side view showing the A-phase 
section on the bottom, with only a few vortices left in the 
outermost ring, while in the B-phase the AB interface at large 
r is covered with a vortex sheet [(|. Above the interface the 
vortices are helically twisted around the central axis. Their 
polarization parallel to the axis increases with distance from 
the interface. Parameters: T = 0.22 T c and P = 29 bar, which 
corresponds to qa ~ 2 and «b = 4.3 x 10 -3 [while a' A ~ 0.8 
and a' B ~ 0], cylinder radius R = 3 mm, the lengths of the 
A and B-phase sections are 5 and 10 mm, the cylinder axis is 
aligned parallel to the rotation axis. 

is not constant across the cross section of the cylinder 
and the superfluid fraction is not in solid-body rotation, 
but relaxes to the final state v s = faster in the center 
than at large radii, (ii) Because of the central vortex-free 
dome, which borders to the AB interface and decreases 
in radius further away at increasing z, the vortices be- 
come curved along the cylinder axis z. Differences in 
the azimuthal velocities as a function of z cause them to 
become helically twisted fioj ]. The inhomogeneous twist 
introduces reconnections and turbulence, which speed up 
the B-phase response. On the AB interface the remain- 
ing vortices curve parallel to the boundary and extend 
radially outward, ending perpendicularly on the cylinder 
wall. This is a B-phase vortex sheet, analogous to the A- 
phase vortex sheet covering the AB interface in rotation 
when the B-phase is free of vortices . 

Spin up, in turn, depends on the number and con- 
figuration of remanent vortices in the B-phase section. 
The influence of the AB interface is to speed up the re- 
moval of remnants. Thus also spin up is different from 
that measured for the single-phase sample in the same 
conditions. These differences in the spin-down/spin- up 
responses of the two-phase and single-phase samples are 
most noticeable at the lowest temperatures. We concen- 
trate here on measurements at 29 bar liquid 3 He pressure 
and 0.20 T c . This is the minimum temperature for the 
current setup, which is limited by a residual heat leak of 
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FIG. 3: NMR absorption spectra during spin down. The line 
shapes (red curves) have been calculated using the cylindrical 
shell model for the two-phase sample and show the counter- 
flow (cf) peak on the right at large frequency shifts / — /l 
from the Larmor frequency /l . The decreasing cf peak heights 
(red lozenges) correspond to increasing radius r a of the cen- 
tral vortex- free cylinder in 20 )j,m steps, starting from r B =0. 
The blue curves are from the measurements where the cf peak 
is monitored continuously by sweeping around its maximum 
(blue dots) during its decay. The green triangles show the 
measured peak trajectory in laminar spin down in the ab- 
sence of the A-phase layer [ill]. Parameters: T = 0.20 T c , 
P = 29 bar, ffc = l.Orad/s, d£l/dt= -0.03rad/s 2 . 

~ 15 pW through the lower orifice of 0.3 mm diameter. 

NMR techniques:-In vortex-free rotation at con- 
stant f2 the B-phase order parameter texture is modified 
by the azimuthally flowing superfluid counterflow (cf) at 
the velocity v = v n — v s = flxr. This solid-body-like cf 
reorients the B-phase anisotropy axis, which is induced 
by the NMR polarization field. The order parameter 
distribution is thereby changed, and a large NMR fre- 
quency shift appears. In the NMR spectrum the shift is 
expressed as a so-called cf peak which can be calibrated 
as a function of the experimental variables f2, T, and P 
[lH . The result is a unique trajectory as a function of 
for the peak absorption vs the frequency shift at fixed T 
and P, which identifies the solid-body velocity distribu- 
tion (Fig. [3]) • During the slow spin down of the superfluid 
fraction after a sudden stop of rotation, when v n = and 
v s arises from the distribution of vorticity, also a cf peak 
is formed. If the spin-down response is laminar, as is 
the case for the single-phase sample, the vortices remain 
approximately straight and maintain a constant density 
across the cross section of the cylinder, i. e. they are in the 
solid-body configuration and the trajectory of this tran- 
sient cf peak matches that measured in the vortex-free 
state at stationary conditions as a function of f2 [1| ■ 

The cf peak recorded during the spin down of the two- 
phase sample follows a different trajectory of peak height 
vs frequency (Fig. [3]). Its frequency shift proves to be 
larger, presumably owing to the compression of vortices 
to large radii around the vortex-free dome in the center. 
This observation suggests a simple model how to analyze 



the changed vortex distribution. Assume that all vor- 
tices are compressed into an outer cylindrical shell with 
the initial vortex density 257o/ K and that the empty inner 
region increases gradually in radius r s . In this model the 
azimuthal velocity of the superfluid component is v s = 
inside the central cylinder r < r s , while in the outer 
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cylindrical shell r B < r < R it is v B 
Here fl(t = 0) = Qq « l.Orad/s is the value extrapolated 
back to t = 0, the moment when the cryostat rotation 
comes to a stop and the cf peak height passes through 
its maximum value. To compare results with and with- 
out the A-phase layer, we use the normalized solid-body 
vortex density equivalent f2/f2o = 1 — (r s /R) 2 . In Fig. [3] 
the NMR spectra have been calculated for every 20 fj,m 
step increase in r s The trajectory of the cf peaks 

calculated in this way agrees closely with that of the mea- 
sured spin-down response. Thus the simple model with 
only azimuthal cf appears to capture the dominant fea- 
tures, although it neglects all dependence on the distance 
from the AB interface and on twisted or tangled vortices. 

Spin-down measurements:-For quantitative analy- 
sis, the decreasing cf peaks at t > were fitted by means 
of the texture calculation procedure to the cylindrical 
shell model with one free parameter, by adjusting the 
inner cylinder radius r s and assuming a vortex density 
outside, 2Qq/k, which is constant as a function of r and 
t. The resulting fit in the inset of Fig. @] reveals a rapidly 
growing central region where the superfluid component 
has already come to rest: in 50 s the radius of the vortex- 
free central region has grown to half of the cylinder radius 
and at this point ~ 25 % of the vorticity has annihilated. 

In Fig.[5]thc spin-down responses of the azimuthal flow 
f2(t)/f2o are compared in the two cases, with and with- 
out the A-phase layer. The measurement is performed 
by decelerating the rotation drive from fij to zero at 
dfl/dt = — 0.03rad/s 2 and the cf peak is recorded as 
a function of time. With decreasing temperature the cf 
peak measurement procedure runs into problems: the cf 
peak height is reduced and reaches zero at higher vortex- 
free flow velocities [ll|. At 0.20 T c in Fig. g] the cf peak 
is lost when the vorticity has dropped by ~ 25 % (filled 
symbols). This explains the blank region where there are 
no data points. In this region more data (open symbols) 
can be retrieved by increasing rotation suddenly back 
to some large value, where the cf peak can be recorded 
while it decays during spin up. This data can then be 
extrapolated back to the moment when rotation was in- 
creased. Combining the two data sets, the spin-down 
response is seen to be a smooth monotonic decay, both 
with and without the A-phase layer. However, in the for- 
mer case it is appreciably faster and not of the laminar 
form fl(t) = Oo/(1 + i/ r -D-)> wncrc a single time constant 
Tjj. = (2aQ )~ 1 fits the data. 

For comparison, Fig. [4] also shows the corresponding 
spin-down response for the calculation in Fig. [5] at dif- 
ferent distances Az from the AB interface. We find that 
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FIG. 4: Normalized spin-down of the azimuthal flow Q(t)/f2o- 
The response of the two-phase sample (blue triangles), ana- 
lyzed as shown in the inset, is compared to that measured in 
the absence of the A phase layer (green circles). Filled symbols 
correspond to cf peaks measured during spin down, while open 
symbols represent extrapolations from subsequent spin ups 
(see text). The broken line curves represent laminar fits with 
fi(t) = Qo/i^ + t/rsj.), where = (2afio) -1 : the dotted curve 
with = 740 s represents the fit to the data measured in the 
absence of the A-phase layer, while the dash-dotted curve is 
a fit to the laminar tail of the two-phase-sample data with 
rj| = 67 s. The solid curves are calculations of the normalized 
azimuthal velocity v s ,4>(R) / (QoR) with the parameter values 
from Fig. [21 The three curves represent different distances Az 
from the AB interface: red, Az — 1 mm; blue, Az = 5 mm; 
green, Az — 10 mm. The dashed fit with = 170 s is the 
laminar tail of the late-time average of the three calculated 
curves. (Inset) Fit to the cylindrical shell model of the two- 
phase sample. (Right vertical axis) Normalized radius r s /R 
(red triangles) of the vortex-free central cylinder, when fitted 
to the measured spin-down as a function of time t. (Left ver- 
tical axis) Equivalent normalized solid-body vortex density 
f2/fio = 1 - (r s /R) 2 (blue dots). 

the result is not very sensitive on Az. Comparing the 
shapes of the response curves we see that in the presence 
of the A-phase layer an abrupt shoulder appears, result- 
ing from faster spin down than dictated by the laminar 
dependence. This is then towards the end of the decay 
followed by a slow final laminar tail. The fast section of 
spin down decay is present in both the calculated and 
measured results in the presence of the A-phase layer. 
This difference from the laminar dependence quantifies 
the turbulent dissipation [l3j |. 

Spin-up measurements:-!!! similar manner we can 
determine from the cf peak response the decay of az- 
imuthal flow during spin up, after a rapid increase of Q, 
Q . We restrict the discussion to the simplest case where 
the acceleration is started from rest (fij = 0) and is fin- 
ished at fif < f^AB (T, P) . In the A-phase section vor- 
tices are formed rapidly and independently during the 
rotation increase. On the AB interfaces they curve radi- 
ally outward, covering the interface as a vortex sheet Q. 
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FIG. 5: Measured cf peak height as a function of time in spin 
up. (Filled symbols) Spin up after having waited for 40 min at 
0.20 T c in zero rotation for remanent vortices from a previous 
spin-down measurement to annihilate. The response is much 
slower in the presence of the A-phase layer (red triangles) 
than in its absence (blue squares). (Open symbols) Spin up 
from a state with ~ 120 remnants. Spin up is now faster with 
the A-phase layer (red triangles) than without (blue squares). 
The vertical arrow denotes the moment when the vortex front 
enters the NMR coil. In these measurements rotation is in- 
creased from zero to flf = l.Orad/s at dfl/dt = 0.03rad/s 2 . 
Time t = marks the moment when fif is reached. The 
number of remnants is estimated from the cf peak height ex- 
trapolated to t = 0. 

Each of these connects to a B-phase vortex across the AB 
interface later, when B-phase vortices are formed. The 
B-phase spin up may proceed in two different ways, de- 
pending on the number and confi gura tion of remanent 
vortices in the initial state at rest [14j. (1) If there is a 
fair number of remnants evenly distributed along the B- 
phase section, then the response might be a slow axially 
homogeneous build up of the vortex density with solid- 
body distribution, which is carried to completion all the 
way to the equilibrium vortex state. In this case the cf 
peak follows the solid-body rotation trajectory in Fig. [3] 
(2) If there are only few remnants, the slow axially ho- 
mogeneous response might be terminated in a localized 
turbulent burst of vortex formation, followed by a subse- 
quent axial expansion of the vortices along the rotating 
cylinder 15( . This spin- up process results in a markedly 
different cf peak response, one that depends on the dis- 
tance of the turbulent burst from the detector coil. 

In Fig. [S] two cases of spin up with and without the 
A-phase layer are compared. Here the remnants are left- 
over vortices from a previous spin-down measurement. 
In one case a waiting period of 40 min is enforced at zero 
rotation, to allow for annihilation. This is calculated 
from the moment when rotation comes to a stop after 
the preceding spin-down measurement. Owing to the 
faster spin-down in the presence of the A-phase layer, 
the number and length of remnants is smaller so that the 
subsequent spin- up is prominently slower ( red filled tri- 
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angles) than in the absence of the interface (blue filled 
squares). In contrast, when the number of remnants is 
adjusted to be roughly equal and relatively large, the 
spin-up response becomes faster in the presence of the 
A-phase layer (red open triangles), owing to a turbulent 
burst of vortex formation at the AB interface which starts 
the axial motion of vortices along the cylinder. The pas- 
sage of the vortex front through the NMR coil is indicated 
by an abrupt decline of the cf peak height [l6| ■ Presum- 
ably the turbulent burst is triggered by the complex re- 
conncction processes which take place in the vicinity of 
the AB interface when A and B phase vortices connect 
across the interface and move at different velocities. 

Spin-down calculations:-Vortex filament calcula- 
tions have proven instructive for analyzing spin-down re- 
sponses, while vortex formation has turned out to be 
problematic in spin- up Qj}j . Our calculation of spin 
down, of which Fig. [5] is a snapshot, makes the following 
two simplifications: all vortices are considered to be simi- 
lar and singly quantized, while the AB interface is treated 
as a plane where mutual friction dissipation changes dis- 
continuously from qa to qb at z — 5 mm. The viscous 
normal component is locked to the rotation drive. At 
low rotation Oi < SIab the calculation describes the dy- 
namics properly, in spite of the simplifications, and re- 
produces the rapid response in the A-phase section, the 
formation of a B-phase vortex sheet on the AB interface, 
and the helically twisted vortices in the B phase, with all 
the currents which follow from these configurations [l7| . 

Fig. |5] summarizes some of the differences in the spin 
down with and without the A-phase section. Most im- 
portant is the much larger frequency of vortex reconnec- 
tions (red curves) and their spatial distribution. Since 
inter-vortex reconnections in the bulk volume feed the 
turbulence, it is instructive to compare them in the two 
cases. Roughly half of all reconnections occurs within 
0.2 mm of the AB interface while the other half is dis- 
tributed relatively evenly as a function of z above the 
AB interface. Radially the reconnections increase steeply 
towards large radii close to the AB interface, while away 
from the interface the increase is slower. In contrast, in 
the absence of the A phase section spin down is lami- 
nar and the reconnections are concentrated to a narrow 
surface layer on the cylindrical wall, leaving the central 
bulk volume reconnection free Thus here we have a 
clear difference - in the former case reconnections occur 
in the bulk volume and on the AB interface, in particu- 
lar, while the surface layer is in the latter case the place 
where the annihilating and reconnecting vortices transfer 
their angular momentum to the walls. 

The comparison on the decay of the total vortex length 
L(t) underlines in a similar manner the differences be- 
tween turbulent and laminar responses (blue curves): 
turbulent spin down in the presence of the A-phase sec- 
tion is non-monotonic with an initial overshoot and a sub- 
sequent rapid decay, which might be characterized with a 
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FIG. 6: Calculated B-phase spin down for the example in 
Fig. [21 emphasizing differences in the response with and with- 
out the A-phase section. (Red) Reconnection rate dN T /dt 
as a function of normalized time Q,ot (solid line). The non- 
monotonic rate peaks sharply at around 70 s and then declines 
rapidly in ~ 400 s to 10 % of the peak value. In compari- 
son, the reconnection rate is orders of magnitude smaller and 
limited to a surface layer in the absence of the A-phase sec- 
tion (dashed line). (Blue) The normalized total vortex length 
L(t)/L(0) also displays an overshoot and peaks at around 
90 s (solid line). It then declines rapidly in 900 s to 10 % of its 
peak value. The later decay is slower with £ _1 dependence. 
In the absence of the A-phase section L(t)/L(0) displays a 
monotonically decreasing curve with slow (1 + t/rjj.) -1 depen- 
dence (dashed line). The parameters are the same as in Fig. [2] 
The reference curves with no A-phase section have been cal- 
culated for a cylinder of 10 mm length, tilted by n — 2° from 
the rotation axis, to break cylindrical symmetry. 

£-3/2 dependence Here part of the kinetic energy of 
the superfluid component is dissipated in reconnections 
and other turbulent excitations. Subsequently the decay 
turns into the slower t^ 1 dependence in the laminar tail. 
In contrast, in the absence of the A phase section L(t) 
decays monotonically with the laminar (1 + t/Tjy) -1 de- 
pendence. Overall we find close correspondence between 
the calculated and measured spin-down responses which 
indicates that differences in vortex velocities across the 
AB interface explain adequately our results. 

Conclusions:-The two-phase superfluid 3 He sample 
provides a unique environment for vortex studies in the 
zero temperature limit, with two interacting supcrfluids 
and their vastly different time scales of vortex flow. This 
leads to unusual vortex configurations in the B-phase sec- 
tion during spin down: a vortex-free dome is formed in 
the center, surrounded by a vortex sheet on the AB in- 
terface and an outer shell of helically twisted vortices 
at high density. Reconnections among the twisted vor- 
tices and on the AB interface concentrate to large radii. 
They bring about increased dissipation and a faster tur- 
bulent spin-down response. These differences from axi- 
ally homogeneous laminar spin down in the absence of 
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the A-phase make it possible to identify the first changes 
in the response from fully laminar towards weakly tur- 
bulent, which arise from the influence of a well-defined 
planar perturbation. The resulting spin-down response 
is different from that when the perturbation is present 
at all surfaces, such as the spin down of a viscous fluid 
in a cylinder |18j, or the spin down of a cube filled with 
superfluid He [5j, |l3| . 
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